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Abstract 

A full quantum mechanical calculation of partial cross-sections leading to different final states 
of antiprotonic helium atom was performed. Calculations were carried out for a wide range of 
antiprotonic helium states and incident (lab) energies of the antiproton. 
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I. INTRODUCTION 



One of the most impressive success stories of the last decade in few-body physics is the 
high precision experimental and theoretical studies of long lived states in antiprotonic helium 
(for an overview see jl|). While the energy levels have been both measured and calculated 
to an extreme precision, allowing even for improvement of numerical values of fundamental 
physical constants, some other relevant properties of these states were studied with consid- 
erably less accuracy. Among these is the formation probability of different metastable states 
(J, v) in the capture reaction 

p + *He^CHe + p) Jv + e-. (1) 

The existing calculations of the capture rates of slow antiprotons in He 0, 0, Q] are based 
on classical or semiclassical approaches and they mainly address the reproduction of the 
overall fraction (3%) of delayed annihilation events. Recent experimental results from the 
ASACUSA project 5|, however, yield some information on individual populations of different 
metastable states, and our aim is to perform a fully quantum mechanical calculation of the 
formation probability of different states in the capture reaction. 



II. CALCULATION METHOD 



The exact solution of the quantum mechanical four-body problem, underlying the re- 
action (^) is far beyond the scope of this work, and probably also of presently available 
calculational possibilities. Still, we want to make a full quantum mechanical, though ap- 
proximate, calculation of the above process. Full is meant in the sense that all degrees of 
freedom are taken explicitly into account, all the wave functions we use, are true four-body 
states. 

The simplest way to realize this idea is to use the plane wave Born approximation which 
amounts to replacing in the transition matrix element the exact scattering wave function 
^ by its initial state $j which preceded the collision. 

Tf i = (9 f \V f \9t)M(* f \Vf\$ i ). (2) 
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In our case the initial and final wave functions were taken in the form: 

$kK,(ri,r 2 ,R) = $ j6 ( ri) r 2 )_I_e iKiR 



f (n. n„ Ti\ = <b, (V e iK /P 



'2 



(2tt) 3 / 2 

where r, are the vectors pointing from helium to the i-th electron, R is the vector between 
He and p, while p { are the Jacobian vectors connecting the electrons with the center of 
mass of the He—p system. For the He the ground state wave function we used the simplest 
effective charge hydrogen-like ansatz [f| 

$ffe(ri,r 2 ) = N exp(-a(n + r 2 )). (3) 

For the antiprotonic helium wave function we used the Born-Oppenheimer form j?l 
which correctly reflects the main features of the final state: 

* JV (p, R) = YMfyMn R) (4) 

where <fii a (p] R) is a ground state two-center wave function, describing the electron motion 
in the field of He and p separated by a fixed distance R, while XJv{R) is the heavy-particle 
relative motion wave function corresponding to ( 4 Hep e~) angular momentum J and "vi- 
brational" quantum number v. The transition potential Vf is obviously 

2 1 1 

V f = + i BT + i 1 

r 2 \Ti - R| |ri - r 2 | _ 

The partial cross-section leading to a certain antiprotonic helium state (Jv) can be written 

as 

^ = W^^^f J d ^f \(*Jv# f \Vf\*He# i )\ 2 (5) 

The angular integrations occurring in the evaluation of Eq. (J3J) were carried out exactly, 
using angular moment algebra, while the 3-fold radial integrals were calculated numerically. 

The general expression (jSJ) for the cross-section leading to a specific state (Jv) can be 
rewritten in terms of matrix element between angular momentum eigenstates as 

a Jv = (2tt) 4 ^ ^ N £(2A + 1) \Mff (6) 

1 A,l 

with 

M% = ( [®j v *k,M)]m I V f I ^(R)]4r >. ( 7 ) 



where (f>K,i{^) denotes free states which definite angular momentum 

•Mr) = J^ji(Kr)Y lm (r) 

and [ ] J M stands for vector coupling. Since the angular momentum of the He ground state is 
zero, the total angular momentum A of the incident side is carried by the antiproton. A given 
antiprotonic helium final state Jv can be formed with different total angular momenta A 
depending on the orbital momentum I carried away by the emitted electron. Our calculations 
show, that only the values I = and / = 1 give a non-negligible contribution to the sum in 
Eq. ©. 



III. RESULTS AND DISCUSSION 



We have calculated the partial population cross-sections aj v for states with angular mo- 
mentum J and energy Ej v in the interval J = 25 — 51, Ej v = —(4.9 — 2.3) a.u. For 
these states all the energetically allowed transitions were calculated for incident antiproton 
energies in the range 5 — 30 eV. 

Our overall results are presented on the contour plots of Fig. [0 The black line separates 
the regions of short-lived and long-lived states. The latter (on the right side of the line) are 
selected according to the usual criterium of Auger-electron orbital momentum I Auger — 4. 

In Figs.|2HHwe tried to illustrate the dependence of certain selected cross-sections o~j v (E) 
on their parameters. In Fig. El we displayed certain cross-sections for various incident en- 
ergies as a function of antiprotonic helium angular momentum J, connecting points, which 
correspond to a certain vibrational quantum number v, while in Fig. El the connected points 
belong to the same principal quantum number N = J + v + 1. 

Fig. 13] shows the dependence of a few cross-sections on the incident energy of antiproton. 

A table, containing all our results, would not fit the size of this paper, how- 
ever, it can be obtained from the authors upon request, or seen/downloaded at 
http : //www.rmki . kf ki .hu/~revai/table .pdf 

It is not easy to draw general conclusions about the relevant physics of antiproton capture 
from this bulk of data. There is, however, a conspicuous feature of our data, which deserve 
some consideration. In a certain region of antiprotonic helium states the dependence of the 
cross-sections on quantum numbers show a smooth, regular pattern, while with increasing 
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FIG. 1: Overall distribution of calculated cross-sections over the quantum numbers J and v for 
different incident antiproton energies. The black line separates the short lived (on the left) and long 
lived states (on the right side of the line) of antiprotonic helium. All cross-sections are measured 
in units of Oq, ao being the atomic length unit. 

excitation energy this behavior becomes irregular. On Fig. ^this is seen as a transition from 
almost parallel stripes to an "archipelago" type structure, while on Figs. I2HH1 the smooth 
lines become oscillatory. 

In order to reveal the origin of this phenomenon we looked into the structure of the matrix 
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FIG. 2: Cross-sections for the lowest few vibrational quantum numbers v and different incident 
antiproton energies. 



element in Eq (j2J) and found that its actual value is basically determined by the integration 
over R. We can rewrite (JJJ) as 



X j v {R)F{R)j A {K i R)RdR, 



(8) 



where XJv{R) is the He — p relative motion wave function in the BO state Jv, j A (KiR) is the 
spherical Bessel function of the incident p, and F(R) contains all rest: the potentials, the 
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FIG. 3: Examples of cross-sections with fixed iV in the "smooth" and "oscillatory" regime (see 
text). 
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FIG. 4: Dependence of selected cross-sections on incident antiproton energy. 



angular integrals, and the integrals over the electron coordinates. The expression (JHJ) can be 
considered as a kind of "radial, one- dimensional" Born approximation for the transition of 
the antiproton from the initial state j\(KiR) to the final state XJv{R) and F(R) plays the 
role of the potential. 

Fig.Elshows a few characteristic plots of j\(KiR), xjv(R), and F(R). It can be seen, that 
F(R) depends very weakly on the quantum numbers of the transition, thus its interpretation 
as "transition potential" is not meaningless. The other essential conclusion from Fig. El is, 
that the value of the integral is basically determined by the overlap of two rapidly varying 
functions, j\(KiR) and XJv(R). While XJv(R) is strongly localized with rapid decay in both 
directions R — ► oo and R — > 0, j\(KiR) is rapidly oscillating for large R and — due to the 
high angular momentum A — strongly decreasing in the direction R — > 0. For increasing 
J, xjv slightly moves outwards, while increasing v (the number of its nodes) makes it more 
and more oscillating. Increasing incident energy moves j\(KiR) inwards. According to 
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FIG. 5: Parts of the integrand of Eq ©: j A (KiR) (solid line), XJv(R) (dashed line) and F(R) 
(dotted line). For details see text. 

these general observations, the "smooth regime" of the dependence of partial cross-sections 
on the energy and quantum numbers, corresponds to the situation, when only the "outer" 
tail of XJv{R) and the "inner" tail of j\(KiR) overlap. For increasing energy (incident 
or excitation) the oscillating parts of xjv and j\ might overlap leading to an irregular 
"unpredictable" dependence on quantum numbers and incident energy. 

This idea can be traced on the last two plots of Fig. |2l The XJv(R) functions for a given 
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v are essentially of the same form, only with increasing J they are pushed outwards into the 
region of oscillations of j\(KiR). For the nodeless v — function this leads to a decrease of 
the integral in Eq. JSJ), while each node of the v ^ functions produces a minimum in the 
cross-section when it penetrates into the region of non- vanishing j\(KiR). 

The three graphs on the right side of Fig. |3J demonstrate this phenomenon on the case 
of the 20 eV v = 1 curve of Fig. It can be seen, how the position of the node of xjv(R) 
relative to the first peak of the j\(KiR) Bessel-function brings about the minimum of the 
cross-section. 

IV. CONCLUSIONS 

To our knowledge, this is the first calculation of the process (0) in which realistic final 
state wave functions were used. Due to this fact we think, that our results concerning the 
relative population of different final states might be reliable in spite of the poor treatment 
of the dynamics of the capture process. As for the absolute values of cross-sections, a more 
realistic dynamical treatment of the reaction is probably inevitable. 

The transition matrix elements are basically determined by the overlap of the BO function 
Xjv{R) and the incident Bessel function j\(KiR) of the antiproton. All the rest can be 
incorporated into a potential-like function F(R), which weakly depends on the quantum 
numbers of the transition. This feature will be probably preserved if a more realistic initial 
state wave function (both for the electrons and the antiproton) will be used. 

The "smooth" regime of the quantum number dependence of the partial cross-section al- 
lows to check the existing two "thumb rules" [l|,B| for the most likely populated antiprotonic 
helium states. One of them states, that the mostly populated levels will have 



while according to the other assumption, the maximum of the capture cross-section occurs 
for zero (or smallest possible) energy of the emitted electron and correspondingly for highest 
excitation energy. From our contour plots of Fig. Q we can conclude, that the maximum 
cross-sections occur along a line, which can be approximated by 




(9) 



v(a,, 



max 



a — b * J 



(10) 
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with a ~ 15 — 20 and b ~ 0.4 — 0.45, depending on incident energy. This observation does 
not seem to confirm any of the " thumb rules" . 

As for comparison of our results with the recently obtained experimental data 5], we 
would like to make two remarks. First, our calculations show, that the cross-sections strongly 
depend on incident antiproton energy. Since the energy distribution of the antiprotons be- 
fore the capture is unknown, the direct comparison with the observed data is impossible. 
Secondly, any observed population data inevitably involve a certain time delay after forma- 
tion and thus the effect of "depopulation" due to collisional quenching. Since this effect is 
absent from our calculation, again, the comparison with experimental data is not obvious. 
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